Studies of nonadiabatic dynamics in the singlet fission processes of pentacene dimer via tensor network method.

Singlet fission (SF) is a very significant photophysical phenomenon and possesses potential applications. In this work, we try to give a rather detailed theoretical investigation of the SF process in the stacked polyacene dimer by combining the high-level quantum chemistry calculations and the quantum dynamics simulations based on the tensor network method. Starting with the construction of the linear vibronic coupling model, we explore the pure electronic dynamics and the vibronic dynamics in the SF processes. The role of vibrational modes in nonadiabatic dynamics is addressed. The results show that the super-exchange mechanism mediated by the charge-transfer state is found in both pure electronic dynamics and the nonadiabatic dynamics. Particularly the vibrational modes with the frequencies resonance with the adiabatic energy gap play very import roles in the SF dynamics. This work not only provides a deep and detailed understanding of the SF process but also verifies the efficiency of the tensor network method with the train structure that can serve as the reference dynamics method to explore the dynamics behaviors of complex systems.

Resonance Enhancement of Vibrational Polariton Chemistry Obtained from the Mixed Quantum-Classical Dynamics Simulations.

We applied a variety of mixed quantum-classical (MQC) approaches to simulate the VSC-influenced reaction rate constant. All of these MQC simulations treat the key vibrational levels associated with the reaction coordinate in the quantum subsystem (as quantum states), whereas all other degrees of freedom (DOFs) are treated inside the classical subsystem. We find that, as long as we have the quantum state descriptions for the vibrational DOFs, one can correctly describe the VSC resonance condition when the cavity frequency matches the bond vibrational frequency. This correct resonance behavior can be obtained regardless of the detailed MQC methods that one uses. The results suggest that the MQC approaches can generate semiquantitative agreement with the exact results for rate constant changes when changing the cavity frequency, the light-matter coupling strength, or the cavity lifetime. The finding of this work suggests that one can use computationally economic MQC approaches to explore the collective coupling scenario when many molecules are collectively coupled to many cavity modes in the future.

Influence of Mode-Specific Excitation on the Nonadiabatic Dynamics of Methyl Nitrate (CH3ONO2).

The impact of mode-specific vibrational excitations on initial-preparation conditions was studied by examining the excited-state population decay rates in the nonadiabatic dynamics of methyl nitrate (CH3ONO2). In particular, exciting a few specific modes by adding a single quantum of energy clearly decelerated the nonadiabatic dynamics population decay rates. The underlying reason for this slower population decay was explained by analyzing the profiles of the excited-state potential energy surfaces in the Franck-Condon regions and the topology of the S1/S0 conical intersection. This study not only provides physical insights into the key mechanisms controlling nonadiabatic dynamics but also shows the possibility of controlling nonadiabatic dynamics via mode-specific vibrational excitations.

Ab Initio Molecular Cavity Quantum Electrodynamics Simulations Using Machine Learning Models.

We present a mixed quantum-classical simulation of polariton dynamics for molecule-cavity hybrid systems. In particular, we treat the coupled electronic-photonic degrees of freedom (DOFs) as the quantum subsystem and the nuclear DOFs as the classical subsystem and use the trajectory surface hopping approach to simulate non-adiabatic dynamics among the polariton states due to the coupled motion of nuclei. We use the accurate nuclear gradient expression derived from the Pauli-Fierz quantum electrodynamics Hamiltonian without making further approximations. The energies, gradients, and derivative couplings of the molecular systems are obtained from the on-the-fly simulations at the level of complete active space self-consistent field (CASSCF), which are used to compute the polariton energies and nuclear gradients. The derivatives of dipoles are also necessary ingredients in the polariton nuclear gradient expression but are often not readily available in electronic structure methods. To address this challenge, we use a machine learning model with the Kernel ridge regression method to construct the dipoles and further obtain their derivatives, at the same level as the CASSCF theory. The cavity loss process is modeled with the Lindblad jump superoperator on the reduced density of the electronic-photonic quantum subsystem. We investigate the azomethane molecule and its photoinduced isomerization dynamics inside the cavity. Our results show the accuracy of the machine-learned dipoles and their usage in simulating polariton dynamics. Our polariton dynamics results also demonstrate the isomerization reaction of azomethane can be effectively tuned by coupling to an optical cavity and by changing the light-matter coupling strength and the cavity loss rate.

Quasi-diabatic propagation scheme for simulating polariton chemistry.

We generalize the quasi-diabatic (QD) propagation scheme to simulate the non-adiabatic polariton dynamics in molecule-cavity hybrid systems. The adiabatic-Fock states, which are the tensor product states of the adiabatic electronic states of the molecule and photon Fock states, are used as the locally well-defined diabatic states for the dynamics propagation. These locally well-defined diabatic states allow using any diabatic quantum dynamics methods for dynamics propagation, and the definition of these states will be updated at every nuclear time step. We use several recently developed non-adiabatic mapping approaches as the diabatic dynamics methods to simulate polariton quantum dynamics in a Shin-Metiu model coupled to an optical cavity. The results obtained from the mapping approaches provide very accurate population dynamics compared to the numerically exact method and outperform the widely used mixed quantum-classical approaches, such as the Ehrenfest dynamics and the fewest switches surface hopping approach. We envision that the generalized QD scheme developed in this work will provide a powerful tool to perform the non-adiabatic polariton simulations by allowing a direct interface between the diabatic dynamics methods and ab initio polariton information.

Molecular design of DBA-type five-membered heterocyclic rings to achieve 200% exciton utilization for electroluminescence.

Achieving high exciton utilization is a long-cherished goal in the development of organic light-emitting diode materials. Herein, a three-step mechanism is proposed to achieve 200% exciton utilization: (i) hot triplet exciton (T2) conversion to singlet S1; (ii) singlet fission from S1 into two T1; (iii) and then a Dexter energy transfer to phosphors. The requirement is that S1 should lie slightly lower than or close to T2 and twice as high as T1 in energy. For this, a scenario is put forward to design a series of donor-bridge-acceptor (DBA) type molecules with 2E(T1) ≤ E(S1) < E(T2), in which the Baird-type aromatic pyrazoline ring is used as a bridge owing to its stabilized T1 (1.30-1.74 eV) and different kinds of donors and acceptors are linked to the bridge for regulating S1 (2.35-3.87 eV) and T2 (2.44-3.96 eV). The ultrafast spectroscopy and sensitization measurements for one compound (TPA-DBPrz) fully confirm the theoretical predictions.

Nuclear gradient expressions for molecular cavity quantum electrodynamics simulations using mixed quantum-classical methods.

We derive a rigorous nuclear gradient for a molecule-cavity hybrid system using the quantum electrodynamics Hamiltonian. We treat the electronic-photonic degrees of freedom (DOFs) as the quantum subsystem and the nuclei as the classical subsystem. Using the adiabatic basis for the electronic DOF and the Fock basis for the photonic DOF and requiring the total energy conservation of this mixed quantum-classical (MQC) system, we derived the rigorous nuclear gradient for the molecule-cavity hybrid system, which is naturally connected to the approximate gradient under the Jaynes-Cummings approximation. The nuclear gradient expression can be readily used in any MQC simulations and will allow one to perform the non-adiabatic on-the-fly simulation of polariton quantum dynamics. The theoretical developments in this work could significantly benefit the polariton quantum dynamics community with a rigorous nuclear gradient of the molecule-cavity hybrid system and have a broad impact on the future non-adiabatic simulations of polariton quantum dynamics.

Insight into Isomeric Effect on the Photoluminescence and Mechanoluminescence of Cyanostilbene Derivatives.

Molecular structures, packings, and intermolecular interactions significantly affect the photophysical properties of organic luminogens. In this work, the photoluminescence (PL) and mechanoluminescence (ML) of two pairs of isomers, 1/2 and 3/4, were systematically explored. The fluorescence of crystals 1c and 4c is much brighter than that of their isomers 2c and 3c, respectively. Only 1c is ML-active among all four molecules. Single-crystal structural analysis revealed that isomerization of a substituent group affected their molecular packing and intermolecular interactions. Stronger intermolecular interaction and intact three-dimensional hydrogen-bonded networks were formed only in crystal 1c, which were essential for preventing slippage of molecular layers and generating ML; the other molecules were either lacking π-π interactions or C-H···π interactions. Theoretical calculation suggested that the energy barrier between the Franck-Condon (FC) structure and minimum energy crossing point (MECP) structure of 2/3 was much lower than that of 1/4. Nonradiative decay channels of molecules 2 and 3 were thus more easily activated, which led to their lower quantum yield.

Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods.

An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.

Heavy-Atom-Free Room-Temperature Phosphorescent Rylene Imide for High-Performing Organic Photovoltaics.

Organic phosphorescence, originating from triplet excitons, has potential for the development of new generation of organic optoelectronic materials. Herein, two heavy-atom-free room-temperature phosphorescent (RTP) electron acceptors with inherent long lifetime triplet exctions are first reported. These two 3D-fully conjugated rigid perylene imide (PDI) multimers, as the best nonfullerene wide-bandgap electron acceptors, exhibit a significantly elevated T1 of ≈2.1 eV with a room-temperature phosphorescent emission (τ = 66 µs) and a minimized singlet-triplet splitting as low as ≈0.13 eV. The huge spatial congestion between adjacent PDI skeleton endows them with significantly modified electronic characteristics of S1 and T1 . This feature, plus with the fully-conjugated rigid molecular configuration, balances the intersystem crossing rate and fluorescence/phosphorescence rates, and therefore, elevating ET1 to ≈2.1 from 1.2 eV for PDI monomer. Meanwhile, the highly delocalized feature enables the triplet charge-transfer excitons at donor-acceptor interface effectively dissociate into free charges, endowing the RTP electron acceptor based organic solar cells (OSCs) with a high internal quantum efficiency of 84% and excellent charge collection capability of 94%. This study introduces an alternative strategy for designing PDI derivatives with high-triplet state-energy and provides revelatory insights into the fundamental electronic characteristics, photophysical mechanism, and photo-to-current generation pathway.

Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations.

The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.

Theoretical Insight Into the Ultralong Room-Temperature Phosphorescence of Nonplanar Aromatic Hydrocarbon.

Purely aromatic hydrocarbon materials with ultralong room-temperature phosphorescence (RTP) were reported recently, but which is universally recognized as unobservable. To reveal the inherent luminescent mechanism, two compounds, i.e., PT with a faint RTP and HD with strong RTP featured by nonplanar geometry, were chosen as a prototype to study their excited-state electronic structures by using quantum mechanics/molecular mechanics (QM/MM) model. It is demonstrated that the nonplanar ethylene brides can offer σ-electron to strengthen spin-orbit coupling (SOC) between singlet and triplet excited states, which can not only promote intersystem crossing (ISC) of S1→Tn to increase the population of triplet excitons, but also accelerate the radiative decay rate of T1→S0, and thus improving RTP. Impressively, the nonradiative decay rate only has a small increase, owing to the synergistic effect between the increase of SOC and the reduction of reorganization energy of T1→S0 caused by the restricted torsional motions of aromatic rings. Therefore, a bright and long-lived RTP was obtained in aromatic hydrocarbon materials with twisted structure. This work provided a new insight into the ultralong RTP in pure organic materials.

Investigation of nonadiabatic dynamics in the photolysis of methyl nitrate (CH3ONO2) by on-the-fly surface hopping simulation.

The photolysis mechanism of methyl nitrate (CH3ONO2) was studied using the on-the-fly surface hopping dynamics at the XMS-CASPT2 level. Several critical geometries, including electronic state minima and conical intersections, were obtained, which play essential roles in the nonadiabatic dynamics of CH3ONO2. The ultrafast nonadiabatic decay dynamics to the ground state were simulated, which gives a proper explanation on the broad and structureless absorption spectra of CH3ONO2. The photodissociation channels, including CH3O + NO2, CH3O + NO + O, and others, as well as their branching ratios, were identified. When the dynamics starts from the lowest two electronic states (S1 and S2), the CH3O + NO2 channel is the dominant photolysis pathway, although we observed the minor contributions of other channels. In contrast, when the trajectories start from the third excited state S3, both CH3O + NO2 and CH3O + NO + O channels become important. Here the CH3O-NO2 bond dissociation takes place first, and then for some trajectories, the N-O bond of the NO2 part breaks successively. The quasi-degeneracy of electronic states may exist in the dissociation limits of both CH3O + NO2 and CH3O + NO + O channels. The current work provides valuable information in the understanding of experimental findings of the wavelength-dependent photolysis mechanism of CH3ONO2.

Multi-Resonance Deep-Red Emitters with Shallow Potential-Energy Surfaces to Surpass Energy-Gap Law*.

Efficient organic emitters in the deep-red region are rare due to the "energy gap law". Herein, multiple boron (B)- and nitrogen (N)-atoms embedded polycyclic heteroaromatics featuring hybridized π-bonding/ non-bonding molecular orbitals are constructed, providing a way to overcome the above luminescent boundary. The introduction of B-phenyl-B and N-phenyl-N structures enhances the electronic coupling of those para-positioned atoms, forming restricted π-bonds on the phenyl-core for delocalized excited states and thus a narrow energy gap. The mutually ortho-positioned B- and N-atoms also induce a multi-resonance effect on the peripheral skeleton for the non-bonding orbitals, creating shallow potential energy surfaces to eliminate the high-frequency vibrational quenching. The corresponding deep-red emitters with peaks at 662 and 692 nm exhibit narrow full-width at half-maximums of 38 nm, high radiative decay rates of ca. 108  s-1 , ≈100 % photo-luminescence quantum yields and record-high maximum external quantum efficiencies of ca. 28 % in a normal planar organic light-emitting diode structure, simultaneously.

Prediction of the excited-state reaction channels in photo-induced processes of nitrofurantoin using first-principle calculations and dynamics simulations.

The understanding of the photochemistry of antibiotic compounds is important because it gives the direct information on the possible environmental pollution caused by them. Due to their large size, the theoretical studies of their excited-state reactions are rather challenging. In current work, we combined the on-the-fly trajectory surface-hopping dynamics, conical-intersection optimizations and excited-state pathway calculations to study the photochemistry of the trans-isomer of nitrofurantoin, a widely-used drug to treat the urinary tract infections. The dynamics-then-pathway approach was taken. First the trajectory surface hopping dynamics at the state-averaged complete-active-space self-consistent-field (SA-CASSCF) level with small active space and small basis sets were run. Second, the minimum-energy conical-intersection optimizations were performed. Finally the excited pathways from the Frank-Condon region to different reaction channels were built at the multi-state multi-reference second-order perturbation (MS-CASPT2) level with large active space and large basis set. Several possible channels responsible for the photo-induced reaction mechanism of the trans-nitrofurantoin were obtained, including the cleavage of the NO bond of the NO2 moiety, the photoisomerization at the central CN bond, and other internal conversion channels. Our findings give some preliminary explanations on available experimental observations. It is also demonstrates that the current theoretical approach is a powerful tool to explore the excited-state reactions in the photochemistry of media-sized or large-sized drug compounds.

Brightening up Circularly Polarized Luminescence of Monosubstituted Polyacetylene by Conformation Control: Mechanism, Switching, and Sensing.

The first example of luminescent monosubstituted polyacetylenes (mono-PAs) is presented, based on a contracted cis-cisoid polyene backbone. It has an excellent circularly polarized luminescence (CPL) performance with a high dissymmetric factor (up to the order of 10-1 ). The luminescence stems from the helical cis-cisoid PA backbone, which is tightly fixed by the strong intramolecular hydrogen bonds, thereby reversing the energy order of excited states and enabling an emissive energy dissipation. CPL switches are facilely achieved by the solvent and temperature through reversible conformational transition. By taking advantages of fast response and high sensitivity, the thin film of mono-PAs could be used as a CPL-based probe for quantitative detection of trifluoroacetic acid with a wider linear dynamic range than those of photoluminescence and circular dichroism. This work opens a new avenue to develop novel smart CPL materials through modulating conformational transition.

On-the-Fly Symmetrical Quasi-Classical Dynamics with Meyer-Miller Mapping Hamiltonian for the Treatment of Nonadiabatic Dynamics at Conical Intersections.

The on-the-fly version of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (SQC/MM) is implemented to study the nonadiabatic dynamics at conical intersections of polyatomic systems. The current on-the-fly implementation of the SQC/MM method is based on the adiabatic representation and the dressed momentum. To include the zero-point energy (ZPE) correction of the electronic mapping variables, we employ both the γ-adjusted and γ-fixed approaches. Nonadiabatic dynamics of the methaniminium cation (CH2NH2+) and azomethane are simulated using the on-the-fly SQC/MM method. For CH2NH2+, both ZPE correction approaches give reasonable and consistent results. However, for azomethane, the γ-adjusted version of the SQC/MM dynamics behaves much better than the γ-fixed version. Further analysis indicates that it is always recommended to use the γ-adjusted SQC/MM dynamics in the on-the-fly simulation of photoinduced dynamics of polyatomic systems, particularly when the excited state is well separated from the ground state in the Franck-Condon region. This work indicates that the on-the-fly SQC/MM method is a powerful simulation protocol to deal with the nonadiabatic dynamics of realistic polyatomic systems.

Analysis of bath motion in MM-SQC dynamics via dimensionality reduction approach: Principal component analysis.

The system-plus-bath model is an important tool to understand the nonadiabatic dynamics of large molecular systems. Understanding the collective motion of a large number of bath modes is essential for revealing their key roles in the overall dynamics. Here, we applied principal component analysis (PCA) to investigate the bath motion in the basis of a large dataset generated from the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian nonadiabatic dynamics for the excited-state energy transfer in the Frenkel-exciton model. The PCA method clearly elucidated that two types of bath modes, which either display strong vibronic coupling or have frequencies close to that of the electronic transition, are important to the nonadiabatic dynamics. These observations were fully consistent with the physical insights. The conclusions were based on the PCA of the trajectory data and did not involve significant pre-defined physical knowledge. The results show that the PCA approach, which is one of the simplest unsupervised machine learning dimensionality reduction methods, is a powerful one for analyzing complicated nonadiabatic dynamics in the condensed phase with many degrees of freedom.

Achieving Pure Green Electroluminescence with CIEy of 0.69 and EQE of 28.2% from an Aza-Fused Multi-Resonance Emitter.

Pure green emitters are essential for realizing an ultrawide color gamut in next-generation displays. Herein, by fusing the difficult-to-access aza-aromatics onto B (boron)-N (nitrogen) skeleton, a hybridized multi-resonance and charge transfer (HMCT) molecule AZA-BN was successfully synthesized through an effective one-shot multiple cyclization method. AZA-BN shows pure green fluorescence with photoluminance quantum yield of 99.7 %. The corresponding green device exhibits a maximum external quantum efficiency and power efficiency of 28.2 % and 121.7 lm W-1 , respectively, with a full width half maximum (FWHM) of merely 30 nm and Commission Internationale de l'Eclairage (CIE) coordinate y of 0.69, representing the purest green bottom-emitting organic light-emitting diode.

Online high-pH reversed-phase liquid chromatography × low-pH reversed-phase liquid chromatography tandem electrospray ionization mass spectrometry combined with pulse elution gradient in the first dimension for the analysis of alkaloids in Macleaya cordata (willd.) R. Br.

An online high-pH reversed-phase liquid chromatography× low-pH reversed-phase liquid chromatography tandem electrospray ionization mass spectrometry combined with pulse elution gradient in the first dimension was constructed to separate and identify alkaloids from Macleaya cordata (willd.) R. Br. The modulation was performed by using a dual second dimensional columns interface combined with a make-up dilution pump, which is responsible for dilution and neutralization of the first dimensional effluent, and the dual second dimensional columns integrated the trapping and the separation function to reduce the second dimension system dead volume. Taking advantage of the dissociable characteristics of alkaloids, mobile phases with different pH values were applied in the first dimension (pH 9.0) and the second dimension (pH 2.6) to improve the orthogonality of two-dimension separation. Besides, the pulse elution gradient in first dimension and second dimensional gradient were carefully optimized and much better separation was achieved compared to the separation with the traditional two-dimensional liquid chromatography approach. Finally, mass measurement was performed for alkaloids in M. cordata (willd.) R. Br. by coupling proposed two-dimensional liquid chromatography system with triple quadrupole mass spectrometry, and 39 alkaloids were successfully identified by comparing the obtained result with the former reported results.